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In this work, we show that reversible photodegradation of Disperse Orange 11 doped in PMMA 
is not due to dye diffusion - a common phenomenon observed in many dye-doped polymers. The 
change in linear absorbance due to photodegradation of the material shows an isobestic point, 
which is consistent with the formation of a quasi-stable damaged species. Spatially-resolved ampli- 
fied spontaneous emission and fluorescence, both related to the population density, are measured by 
scanning the pump beam over a burn mark. A numerical model of the time evolution of the popu- 
lation density due to diffusion is inconsistent with the experimental data suggesting that diffusion 
is not responsible. 



I. INTRODUCTION 



Photodegradation in dye doped polymers was exten- 
sively studied and modeled in several systems [l|, 
Recently, reversible photodegradation was observed for 
dye molecules when doped in polymers @ but not in 
hquidsQ- A logical hypothesis is that dye diffusion is 
responsible for this effect. 

Dye diffusion is a common phenomenon in azo dye- 
doped polymer systems [3]. There are many applications 
that are based on dye diffusion such as thermal trans- 
fer printing [sj or holographic surface relief gratings [6]. 
Photoinduced dye diffusion is observed in various dye- 
polymer systems and has been thoroughly studied over 
the past decade 

We previously reported on candidate mechanisms that 
could be responsible for reversible photodegradationQ 
in l-amino-2-methyl anthraquinone (commonly known 
as Disperse orange 11) when doped in Poly (methyl 
methacry late) (PMMA). This paper focuses specifically 
on eliminating the likelihood of dye diffusion alone as a 
cause of self-healing effects. Anderson et aVs imaging 
studies of photodamage show the visual recovery of a 
burn track, yet it alone does not fully rule out the possi- 
bility of diffusion. 

In this paper, we present experimental observations of 
changes in absorption spectrum of a damaged area, which 
support the hypothesis of the formation of a damaged 
species and subsequent self healing back to the original 
molecule. We also report on measurements of the time 
evolution of the population density profile and show that 
they are inconsistent with theoretically-simulated profiles 
due to diffusion. 
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FIG. 1: Isobestic behavior in linear absorbance suggests 
molecular structure changes. Sample used: 6g/l of disperse 
oragnge 11 doped in PMMA made into a thin film of thick- 
ness of the order of lOO/xm. Inset shows a cartoon of pump 
line (dark orange) and probe (white circle) overlapping in the 
sample. 



II. EXPERIMENTS 

A frequency doubled Nd:YAG nanosecond pulsed laser 
at a wavelength of 532 nm is brought to a line focus 
of 100/Ltm thick on a sample to induce optical damage. 
This pump beam also produces amplified spontaneous 
emission (ASE) light that propagates along the excitation 
line. The absorption spectrum is measured as a function 
of time using a white light probe that spatially overlaps 
the pump beam at the sample. The inset in Figure [T] 
shows a cartoon of the overlap between the probe beam 
and pump beam line. 

The sample is pumped for 20 minutes, during which 
time the absorption spectrum is measured at t = 0, 10 
and 20 to probe photodegradation. Subsequently, the 
laser is blocked to allow the sample to self heal. The 
absorption spectrum is measured again after 2 hours of 
healing. The experiment is described in more detail in 
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FIG. 2: Fluorescence and ASE intensity as a function of ex- 
citation line position. 



the literature 

Diffusion is characterized by a decrease of population 
in the bright regions, leading to an absorption spectrum 
of the same shape but decreased area. Figure 1 shows the 
absorption spectrum plotted as a function of energy at 
various times with the initial spectrum subtracted. The 
dotted curve shows the initial absorption spectrum of the 
sample. 

An isobestic point for a series of linear absorption spec- 
tra is an indication of the conversion of one molecular 
structure into another one, or changes in the molecular 
structure; but, is not expected for diffusion. The de- 
pletion of undamaged molecules is clear from the dip at 
2.57eV that characterizes the DOll molecule, while the 
growing peak at 3.2eV is an indication of the formation 
of a damaged population. There is a clear isobestic point 
between these two regions. After two hours, the damaged 
population characterized by the peak at 3.2eV decreases 
and the dip at 2.57eV decreases showing recovery of the 
sample. A two peak fit to the curves (not shown) ap- 
proximately yields an equal area under the dip and peak 
showing the conservation of population of molecules in 
the region probed with the white light source. While this 
data clearly shows the conversion between two molecular 
structures, it does not rule out contributions from diffu- 
sion. 

The amplified spontaneous emission (ASE) intensity 
and the fiuorescence intensity are related to the popula- 
tion density in the pump region. As such, a mapping of 
the collected fluorescence or ASE light as a function of 
pump position can be used to determine the population 
profile as a function time. 

We use the line focus of the pump laser both to burn a 
damage line into the sample and to measure the popula- 
tion profile by monitoring fiuorescence/ ASE while trans- 
lating the pump line through the burn line. Since the line 
focus is of transverse gaussian intensity profile, we expect 



FIG. 3: ASE scan over a burn mark as the sample recovers. 
Top: the experimental data at several times, and bottom: 
gaussian fits to the experimental data. 



the burn mark to also have a gaussian profile. The fiu- 
orescence/ ASE profile will thus be a convolution of the 
undamaged population density profile and the pump in- 
tensity profile. If diffusion is a dominant mechanism, 
we would expect a buildup of population just outside 
the burn line due to photothermally-induced transport 
of molecules away from the hot area. Figure [5] shows 
that there is no significant increase of fluorescence/ ASE 
emission when the sample is pumped near the edges of 
the burn line. 

Since ASE is a nonlinear process, the width of the cross 
correlation function determined from a scan is narrower 
than it is for linear fluorescence - as observed in Figure 
[21 Figure [3] shows the time evolution of the ASE profile 
as a function of time during recovery. The decrease in 
the width of these curves is consistent with the fact that 
ASE is a nonlinear process. 



III. THEORETICAL MODELING 

In this section, we present a model of dye diffusion in 
a polymer due to a photo-induced temperature gradient. 
We consider the contribution of diffusion to the mecha- 
nisms of photo-degradation and self-healing by simulat- 
ing both the population loss in the illuminated region 
and the return of population when the sample recovers 
in the dark. 



A. Diffusion from a temperature gradient 

It is well known that the equilibrium population den- 
sity of a molecular species is altered by a tempera- 
ture gradient due to molecular diffusion in the host 
material [ill ll^ . Fick's first law for particles in a ma- 
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terial with a temperature gradient is 



J = CMF (r, C)-D (r, C) VC, 



(1) 



where J is the diffusion flux, C is the concentration, M 
is the mobility, F is the force from the temperature gra- 
dient, D is the diffusivity, t is the time, and r is the 
position. 

The force on dye molecules from the temperature gra- 
dient, F, can be expressed in terms of the temperature, 
T, and the transport heat, Q. This gives 



F = -— VT. 
T 



(2) 



The mobility of the dye molecules in the polymer is given 

by 



(3) 
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FIG. 4: The temperature profile from a pump beam that is 
infinitely long in the j/-direction and has a gaussian intensity 
profile in the a;-direction. 



where fi is the chemical potential. Note that this assumes 
a constant volume. We now can rewrite Fick's first law 
for constant volume with concentration and temperature 
gradient as 



J = - ( CMj;VT + D (r, C) VC 



(4) 



The parameters used in each simulation are chosen ac- 
cording to the experimental conditions. The samples are 
approximately 100 /im thick and 1 cm square. The pump 
beam, of gaussian profile, is focused by a cylindrical lens 
onto the flat surface of the sample. The line focus spans 
the entire length of the sample in one direction. The 
transverse beam intensity profile is much smaller than 
the width of the sample. Since measurements are taken 
in the center of the sample, we approximate the beam as 
an infinitely long heat source. 

Next, we evaluate Fick's second law for a constant 
temperature gradient. The constant temperature gra- 
dient approximation due to a photo-induced heat source 
is valid for any system that exhibits a thermal diffusivity 
that is much greater than the diffusivity of the molecules 
through the medium, which is valid for a dye-doped poly- 
mer [ij, [iJl . The diffusion equation for a material with 
a known temperature profile, T{x, y), is 



'dt 



Q 

CM—VT 
T 



D (r, C) VC 



(5) 



A uniform temperature distribution in the y-direction 
from an infinitely long heat source gives 



dy 



(6) 



Figure 2] shows the temperature profile throughout the 
material for an ambient temperature of 297 K and a peak 
temperature of 307 K, values that are typical for the 



pump powers used. Thus, for an initial uniform con- 
centration 



dy 



0. 



(7) 



Therefore, the concentration only varies in the x- 
direction for a line focus heat source along y. 
The diffusion equation in one dimension is 



dC d /^^OdT ^^dC 



(8) 



Assuming that the diffusivity is constant with respect 
to the small changes in concentration. Equation [5] can 
be rewritten using the quotient rule in terms of separate 
derivatives of the concentration and temperature 
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(9) 



where D, M, and Q are approximated as constants with 
respect to the concentration and the spatial coordinate. 
The temperature profile is 



T = To + 



(10) 



where Figure 2] shows the temperature distribution in the 
plane of the thin-film sample. 

Substituting Equation ITOl into Equation HI and setting 



T' — T 



dC 
'dt 



= D 



/To, gives 



dx^ 



MQ 



(11) 
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To avoid small errors at the boundaries due to numerical 
approximations, we impose Dirichlet (rather than Neu- 
mann) boundary conditions at a: = ±L/2, where L is 
the length of each edge of the square sample and the 
boundaries are far from the heat source. This is espe- 
cially important for computing small changes over long 
periods of time. The Dirichlet boundary conditions are 



(12) 



We also assume that the concentration is uniform before 
the laser is turned on, or, 



C (x, t = 0) Co 



(13) 



We begin by first solving the diffusion equation while 
the pump laser is on to determine the initial concentra- 
tion profile when the laser is turned off and the sample 
begins to heal. At this time. Equation [TT] reduces to the 
simplest form of the diffusion equation, 



dC_ 
'dt 



D 



(14) 



Equation [TT] has no transient solutions in an analytical 
form of common functions. Therefore, we must find a 
numerical approximation to determine the concentration 
as a function of time and position due to the tempera- 
ture gradient induced by the laser. The Crank-Nicolson 
method is used to solve all diffusion equations. 

We investigate two cases. The first involves pure dif- 
fusion, which assumes that there are no chemical reac- 
tions. The second assumes that there are two process; 
namely, a partial burn due to photochemistry and dye 
diffusion. For both cases, we use a nonlinear fitting 
scheme to approximately determine the parameters MQ, 
Z), and a in the diffusion equation. Note that we assume 
T' — 10/297 K - values that are typical in our experi- 
ments. This is a good approximation because 



T' 



1 ^ Ji/g-a;2/2cr2 



(15) 



when T' is small, which is the case in our experiments. 

Previous imaging experiments by Anderson et aZ[lo| 
studied the transmission of light in a thin sample of 
DOll-doped PMMA as a function of time while the ma- 
terial is recovering. Assuming that the recovery process 
is due to diffusion, these images can be used to determine 
the fractional change in the concentration as a function 
of position and time. To do so, the Beer-Lambert law is 
inverted to give 



— = 1-— ln(^^^ 

Co abso I Tb 



(16) 



where Co is the initial concentration and abso is the ab- 
sorption coefficient at the peak wavelength of the light 
source used to construct the image. For a 100 /im thick 
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FIG. 5: The concentration fraction from the gaussian fits 
of experimental data(points) plotted along with the diffu- 
sion results(curves). We find that a — 24.4 /xm"^, D = 
0.767 ^im^min"\ MQ = 23.6 ^im^min"\ 



sample at an initial concentration of 9 g/L and constant 
absorption cross section, abso = 2.56. T is the transmit- 
tance of the sample, and Tb is the baseline transmittance 
that is measured prior to illumination with the pump 
light that is used to damage the material. 



B. Parameters for pure diffusion 

The Gauss-Newton method was used to fit the out- 
put data generated by the Crank-Nicolson method to the 
concentration fraction determined from the experimental 
data. The elements of the Jacobian were estimated for 
the iterative fit of the partial differential equation. 

Figure [H] shows a plot of the measured concentration 
fraction profile (symbols) and the concentration fraction 
calculated from Equation [TT] (curves) as a function of 
time using as the initial condition the image of the dam- 
aged sample at the start of self healing. For pure diffu- 
sion, the concentration distribution broadens over time. 
However, the data clearly shows that only the amplitude 
changes, not the width. Thus, we conclude that a pure 
diffusion mechanism cannot explain the data. 



C. Parameters for diffusion and photochemistry 

We now consider the case where there are two pro- 
cesses: a permanent photochemical burn and dye diffu- 
sion. The permanent burn is identified by the residual 
change in transmittance at 1500 minutes. The image is 
found to be constant for up to 4500 minutes. Thus, we 
subtract this residual transmittance to isolate what we 
hypothesize is the diffusion component, which we ana- 
lyze using Equation 1161 
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FIG. 6: The concentration fraction from the gaussian fits 
of experimental data plotted against the diffusion results. 



We find that a 
18.5 ^m^min"^ . 



24.0 /im~\ D = 5.37^im^min"\ MQ 



The new base-lined data is shown in Figure [51 which 
is fit to the diffusion model. If the sample were to relax 
due to diffusion after the laser is turned off, the diffusion 



equation predicts that the width of the concentration pro- 
file increases. Again, we find that the data and theory 
disagree, so that diffusion does not appear to contribute. 



IV. CONCLUSIONS 

We have considered two possible scenarios: pure diffu- 
sion and diffusion with permanent chemical changes due 
to photodegradation. In both cases, the theory of dif- 
fusion and the experimental results do not agree. Thus, 
the contribution of diffusion to self healing of DO 11 in 
PMMA polymer must be negligible. 

The isobestic point and the evolution of the linear ab- 
sorbance spectrum of DO 11 in PMMA suggests that the 
damaged species is structurally distinct from the initial 
molecule. While diffusion due to localized laser heating 
may be present, our measurements show that its effect 
is negligible. Ruling out dye diffusion will allow other 
mechanisms to be tested in the quest towards developing 
an understanding of the self-healing phenomena. 
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